Interplay between phase defects and spin polarization in the specific heat of 
the spin density wave compound (TMTTF) 2 Br in a magnetic field 
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Equilibrium heat relaxation experiments provide evidence that the ground state of the commensu- 
rate spin density wave (SDW) compound (TMTTF)2Br after the application of a sufficient magnetic 
field is different from the conventional ground state. The experiments are interpreted on the basis 
of the local model of strong pinning as the deconfinement of soliton-antisoliton pairs triggered by 
the Zeeman coupling to spin degrees of freedom, resulting in a magnetic field induced density wave 
glass for the spin carrying phase configuration. 

PACS numbers: 75.30.Fv,75.40.Cx,75.10.Nr 



Metastability results from energy minima in a system 
phase space, separated from each other by energy barri- 
ers. In some complex systems such as spin glasses, the 
number of metastable states and the scaling of the energy 
barriers with the system size are such that ergodicity is 
broken, meaning that time averaging is not equivalent to 
ensemble averaging because the system is "trapped" in 
a valley of the energy landscape. Metastability is also 
found in model systems such as molecular magnets or 
Josephson junctions, described by an energy potential 
V(ip) as a function of a single degree of freedom ip. We 
investigate below experimentally and theoretically the 
residual degrees of freedom of a spin density wave (SDW) 
at very low temperature in a magnetic field, interpreted 
as the properties of a classical potential V(<p(yi)) for the 
SDW phase <p(yi) at the coordinate along the chain of 
a strong pinning impurity. 

Below the Peierls transition temperature, charge den- 
sity waves (CDWs) and SDWs in quasi-one-dimensional 
(quasi-lD) compounds are characterized by a spatial 
modulation of the electronic density (or spin) along the 
chains. The phase profile is the result of a compromise 
between the elastic energy that penalizes large phase gra- 
dients, and the pinning energy that tends to fix the phase 
at the strong pinning centers. These two ingredients of 
the Fukuyama-Lee-Rice model lead to metastability 
as the result of collective pinning. Collective pinning is 
however frozen below a glass transition of order ~ 50 K, 
as shown by dielectric susceptibility experiments 0. At 
very low temperature, the residual degrees of freedom in 
a zero magnetic field correspond to the local defects of 
the local model of strong pinning 0, 0> IE IS 13 • Larkin 
and Ovchinnikov have shown that a single strong 
pinning impurity leads to a bound state of an electron- 
like soliton and a hole-like antisoliton. This results in a 
potential V(ip(yi)) with multiple minima 0, leading to 
slow relaxation in agreement with the very low tempera- 
ture heat relaxation experiments 0. 



In order to explore the role of a magnetic field, we 
choose the compound (TMTTF^Br with a sufficiently 
narrow spectrum of relaxation times because of its com- 
mensurate (antiferromagnetic) ground state llflj . This 
allows a systematic study of the equilibrium energy re- 
laxation over almost one decade in temperature and a 
direct comparison to the local model of strong pinning 
without introducing an additional time scale related to 
ageing 0[ll|]. 

The specific heat of about 60 mg of a (TMTTF) 2 Br salt 
was measured at CRTBT-Grenoble in a dilution cryostat 
set-up under magnetic field up to 7 T. This compound 
was previously investigated in zero field in a similar tem- 
perature range . We use a standard relaxation method 
but in contrast to previous experiments, the spe- 
cific heat is determined at equilibrium, once the heat re- 
laxation regime AT(t) as a function of time t, obtained 
after a long heat input, becomes exponential. This in- 
ternal relaxation time t{ n reaches up to more than 10 4 s 
at the temperature T = 60 mK and at high field. This 
technique requires probing times of about 3-^4i in for the 
determination of the time constant r eq of the exponential 
relaxation for t > £j n , and hence extremely stable field 
supply and temperature regulation. The equilibrium spe- 
cific heat C eq reported on Figs. Q and El is obtained from 
the relation C eq = T eq /Rhi, where Rhi is the heat link re- 
sistance between the sample and the regulated cold sink, 
measured at equilibrium under permanent power supply 
(more details will be published separately). 

The temperature dependence of the specific heat C P (T) 
is shown on Fig. ^ for different values of the magnetic 
field. We find experimentally a specific heat propor- 
tional to 1/T 2 , with a prefactor varying by almost two or- 
ders of magnitude when the magnetic field increases from 
/i = 0Tto/i = 7T. This temperature variation is inter- 
preted as the high temperature tail of a Schottky 
due to two level-like systems. The specific heat of a con- 
centration A of two-level systems with an energy splitting 
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FIG. 1: (Color online.) Temperature dependence of the spe- 
cific heat for various magnetic fields. (□, purple) corresponds 
to h = T in a zero field cooled sample, and (V, red) corre- 
sponds to h — T in a zero field cooled sample where the mag- 
netic field was cycled according to ft = T ^ 5.5 T -> T. 
The solid lines are a fit to C p ° (/jft/feT). 



A is C { p ] (A/k B T) = A(A/k B T) 2 /[2cosh(A/2k B T)} 2 , 
that behaves like c£ 0) (A//fc B T) ~ A(A/2k B T) 2 for k B T 
large compared to A. The temperature dependence of 
the specific heat is well described by Cp°\nh/k B T), with 
a splitting A = fih. The value of /x ~ 0.011/iB (with 
/is is the Bohr magneton of an electron) deduced from 
the fit on Fig. ^ is approximately one order of magni- 
tude smaller than the SDW amplitude measured by NMR 
[lij l . This is because solitonic excitations involve a distri- 
bution of spins with a staggered orientation, with a net 
magnetic moment smaller than the moment of an indi- 
vidual spin (see Fig. [3J|. We have chosen the fit leading 
to the smallest value of A ~ 60 J/mol K compatible with 
experiments. This results in a huge number of defects 
induced by the magnetic field (approximately 8 defects 
per unit cell) that cannot be explained by impurities . 
Here we relate the magnetic field specific heat to spin de- 
grees of freedom in a density wave (DW) glass, not to 
the phase excitations of bisolitons in the local model of 
strong pinning where the number of phase defect would 
be equal to the number of pinning centers. The fit of 
the temperature dependence of the specific heat in zero 
field leads to a concentration of defects of approximately 
20%. This concentration is consistent with the Ovchin- 
nikov estimate in (TMTSF) 2 PF 6 @, but is however too 
large to be ascribed to extrinsic impurities. 

We report for the first time metastability induced by 
the magnetic field at a fixed temperature T = 92 mK 
(see Fig- - The sample that has not "seen" the mag- 
netic field follows the branch Cp(h) and remains on 
branch I if the applied magnetic field does not exceed 
h c ~ 5.5 T. If the magnetic field increases above h c , the 
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FIG. 2: (Color online.) Magnetic field dependence of the spe- 
cific heat at a fixed temperature T = 92 mK. The magnetic 
field is first cycled from ft = 0Ttoft = 7T, starting from 
a zero field cooled sample (•, green, branch f), and next de- 
creased back from ft = 7Ttoft = 0T ■ red, branch 2). 
The system remains on branch 2 after additional magnetic 
field cycles. Recovery branch 1 starting from branch 2 (a) 
is explained in the text . The solid (red) line "1" is a fit 
to Cp (nh/ksT) (parameters as on Fig. 0. Branch 2 satu- 
rates to Cg for h <C h c . The spin- up bisoliton contribution 
(dotted lines) vanishes at h ~ h c . 



(2) 

specific heat follows Cp '(h) upon decreasing the mag- 
netic field below h c , and remains on branch 2 if the 
magnetic field is further cycled above h c . This signals 
that a magnetic field h > h c induces a new ground state 
that we identify below as a density wave glass. Branch 
I is recovered with the following history: i) the sam- 
ple is reheated in zero field from 90 mK up to above 
20 K ii) once the sample is cooled down, the magnetic 
field is increased up to 1.5 T and decreased again to 
0.2 T. The raw experimental data are well fitted by 
C { P h2) (h,T) = \c { h2) +C P [a) (nh/k B T)\ mJ/mol K for 
branches 1 and 2, where h is in Tesla and T in Kelvin, 



and with = 120 mJ/mol K, and C^' = 404 mJ/mol 



(2) 



K. We find C^ 0) (nh/k B T) 
k B T much larger than jih. 



190(h/T) 2 mJ/mol K for 



The internal equilibrium time ti n follows an activated 
regime t in — r exp (Ea/T) with an activation energy 
Ea — 0.50 K and a magnetic field-dependent attempt 
time To of the order of a few seconds (see Fig. |3J). In 
zero field, we find an excellent agreement with a previous 
determination of the activation energy, obtained from the 
spectrum of relaxation times . The activated behavior 
is in agreement with the local model of strong pinning 
that we consider now. 

The effective ID Hamiltonian for the phase ip(y) of the 
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FIG. 3: (Color online.) Ahrrenius plot tin versus 1/T of the 
internal relaxation time. The symbols correspond to the same 
values of the magnetic field as on Fig. The experimental 
data are fitted by ti n ~ To(h) exp (Ea/T), with Ea = 0.5 K 
and r (2T) = 4.5 s, r (3T) = 6 s, r (4T) = 10 s, r (5.5T) = 
18 s, r (7T) = 27 s, r (0T) = 0.7 s (branch I), r (0T) = 2.8 s 
(branch 2). 



SDW takes the form fl H S H IS S H : 
H=^Jdy(?^) +wjdy[l-co B <p(y)], (1) 

where y is the coordinate along the chain, vp the Fermi 
velocity, w the commensuration potential. The excita- 
tions in the absence of impurities are pairs of ±2-77 soli- 
tons and antisolitons, the phase of which winds by ±2n 
within a correlation length £ = \Jhvp /2nw. A SDW is 
viewed as the superposition of two out-of-phase CDWs 
for spin-up and spin-down electrons, and the SDW pin- 
ning energy is obtained to second order as fl~l. Hil liH 

n 



imp 



J2i Vi cos [2 (Qyi + ip(yi))] 1 where the sum runs 
over all strong pinning impurities at position y{ along 
the chain. The phase field is ip(y) = f-\{y) — ¥>j,(j/) + t, 
where <p\(y) and <fi(y) are the spin-up and spin-down 
phase fields. The residual specific heat at low field of 
120 mJ/mol K in Cp X \h,T) is related to metastability 
due to spinless bisolitons (bound state of a soliton and 
an antisoliton) |j. The phase profile fb{y) of a bisoli- 
ton is |4| 



u „i's?M) =tal tt |( , sll 



\y-yi\ 



(2) 



The ground state of a bisoliton in the absence of a mag- 
netic field is at energy Eq, separated by a barrier from a 
metastable state at energy E + AE, with AE — AEs = 
16u>£, where Es is the energy of a ±27r soliton in the 
ure system. This defines an effective two-level system 



A Zeeman coupling to the magnetic field is included 
now through Hh = —\xh j dy [pi(y) — Pi(y)], where the 
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FIG. 4: (Color online.) (a): A spin polarized bisoliton made 
of the superposition of a spin-up electron-like soliton in the 
field <pt(v) an d a spin-down hole-like antisoliton in the field 
'Pi in) ' s deconfining (according to the arrows). The dashed 
lines are the continuation of the two solutions, without physi- 
cal meaning, (b): The solitons and the antisolitons are decon- 
fined, and a new bisoliton is generated. In the spin polarized 
case, the energy of (b) is equal to to the energy of (a) minus 
the spin energy 4nhpoQ~ . The impurity is indicated by the 
black box at position yt along the chain. 



electronic density p a (y), with a =t, 4, is defined by 



pa{y) = P0 



dy 



-hpicos (Qy + <Pa{y)), (3) 



where Q is the SDW wave- vector, p is the electronic 
charge per unit length in the absence of deformation. 
The term containing p\ is relevant to local pinning, but 
averages to zero in the Zeeman energy integrated over y 
because of the short scale oscillations at the Fermi wave- 
length. 

A 27r soliton in the spin-up field <p-\(y) and a — 27r 
soliton in the spin-down field <Pi(y) both carry a net 
spin-up because their energy decreases by 2irphp^Q^ 1 
in a magnetic field h. In the absence of impurities, the 
SDW ground state is unstable against the nucleation of 
pairs of solitons and antisolitons carrying a net spin- 
up if the gain in the Zeeman energy exceeds the soli- 
ton energy Es, a condition equivalent to h > h c , with 
h c = 2QEs /ttppq. An upper bound to the temperature 
of the maximum of the Schottky anomaly in a zero field 
is T max ~ 10^20 mK, leading to AE S ~ 2.4T max , so that 
h c — 4Es/p is lower than ~ 2.2-^4.3 T is compatible with 
the cross-over magnetic field obtained experimentally. 

Now, we show how deconfined pairs of solitons and 
antisolitons are nucleated from the strong pinning impu- 
rities. A spinless bisoliton corresponds to <£]{y) = fb(y) 
and 7T— (fi(y) = for all values of y (see Eq. (J2J for Pb{y))- 
Alternatively, the field <pi(y) can be excited: <fi(y) = 
and 7T — <Pi(y) = fb{y) for all values of y. On the other 
hand, a maximally spin polarized bisoliton is obtained 
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FIG. 5: (Color online.) Energy potential (in K) as a function 
of the phase at the position yi of a strong pinning impurity, 
(a): for h = 0, and (b): for ft = 7.3 T (or fih = 0.08 K). 
The pinning energy is such as to reproduce the experimental 
activation energy Ea — 0.5 K (see Fig. |jy. The figures show 
the potential energies of a bisoliton with no additional soli- 
ton/antisoliton pair (solid line, red), with a soliton for y < yi 
and an antisoliton for y > yt (short dashed, blue), and with 
an antisoliton for y < yi and a soliton for y > yi (long dashed, 
green) . 

by reversing the spin in the part of the soliton corre- 
sponding to y > yi, in such a way that <p-\{y) — <Pb{y) 
and 7T — fi{y) — for y < yi, and ip-\(y) — and 
ir — ipi (y) = ifib(y) for y > yi- A sequence generating a de- 
confined soliton-antisoliton pair starts from the minimum 
(1) with a small value ofipi/2ir (see Fig.EJ. The phase ipi 
at the position of the impurity crossed over to the other 
minimum (2) with ipi/2ir ~ 1, and the soliton-antisoliton 
pair deconfines by relaxing to the minimum (3). The 
process can then be iterated, inducing the same random 
configuration of the SDW phase as in a DW glass [3j . The 
spin degrees of freedom on the solitons and antisolitons 
define two-level systems with a splitting proportional to 
the applied magnetic field, which explains the high-field 
contribution Cp°\fj,h/kBT) proportional to (ph/hsT) 2 
(see Fig.|2J). 

The bisoliton contribution in low field (see Fig. [5J is 
obtained by assuming a population of spinless and spin 
polarized bisolitons in thermal equilibrium. The spe- 
cific heat of the spin-up bisolitons with a magnetic mo- 
ment parallel to the applied magnetic field vanishes at 
the magnetic field h c (see Fig- EJ- For this value of the 
magnetic field, the magnetic energy is opposite to the 
Larkin-Ovchinnikov level splitting in zero field, resulting 
in a degenerate two-level system with a vanishingly small 
specific heat. The spinless bisolitons restore a finite spe- 
cific heat (see the total contribution of the bisolitons on 
Fig. 0). If the field is reduced from a large value, the 
pairs of solitons and antisolitons annhilate reversibly for 
h > h c , and the phase is trapped at h = h c . The random 
phase pattern induced by the magnetic field persists if the 
field decreases back to zero. The spin degrees of freedom 
on the residual phase defects explain the specific heat of 



branch 2, larger than the specific heat of branch 1. 

To conclude, we have found experimentally a very per- 
sistent metastable branch in the magnetic field depen- 
dence of the specific heat with a larger specific heat than 
in the zero field case. We interpreted this observation in 
terms of the local model of strong pinning coupled to a 
Zeeman field, in which we find an instability of the SDW 
ground state with a flat phase towards a DW glass with 
a random phase configuration. The additional contribu- 
tion to the specific heat when coming back to zero field 
is explained by the spin entropy due to the magnetic 
moments on the phase defects of the DW glass. Inter- 
estingly, a moderate pressure induces a transition to an 
incommensurate phase in the same compound ; which 
was explained by the formation of discommensurations. 
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